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Abstract. The broad line profiles (H /?) from a nonplanar 
accretion disk resulting from reprocessing of hard X-ray 
radiation are studied in this paper. The simplest model of 
a disk with small inclination and a central point-like source 
of ionizing radiation is adopted. The profiles obtained de- 
pend significantly on the line of sight - they are asymmet- 
ric as well as frequency-shifted in most cases. Since such 
frequency shifted asymmetric profiles are generally not ob- 
served in AGN, the most plausible interpretation may be 
that either the major part of the broad emission does not 
originate from the disk, but from different matter near the 
hole, or the disk is in general planar. The latter situation 
might apply to non-rotating black holes, for which the ac- 
cretion disk is not warped. The emission lines are then 
similar to those from planar disks, i. e. almost symmetric, 
double-peaked and non frequency-shifted. Symmetric and 
double-peaked profiles might also be observed in case of an 
accretion disk transparent to visual light and/or X-rays. 
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1. Introduction 

It is widely believed that Active Galactic Nuclei (AGN) 
are powered by accretion onto a supermassive black hole 
(Rees 1984, Blandford 1990, Ho 1998). This accretion 
probably operates through a geometrically thin disk, 
which is necessary to ensure enough radiation efficiency 
of the process (Rees 1984). Viscous release of energy in 
the disk, due to differential Keplerian rotation, is the 
most commonly invoked mechanism to explain the ob- 
served AGN luminosities, which are usually about IQ''^""''^ 
ergs~^. Through this mechanism, which generates mainly 
visual, UV, and soft X-ray continuum, up to 0.30 of the 
total mass-energy of the flow can be extracted (Thorne 
1974). The observed hard X-ray radiation is probably pro- 
duced by a compact source, located near the center (at 
i? = 10 - 20i?G, where Rq = GMBHC"^ A/bh is the 
black hole mass). The compactness of this X-ray source 
follows from its rapid variability in many objects. 



A thin accretion disk in AGN could not only be the 
energy source but also act as a fluorescent screen, produc- 
ing a significant part of the broad emission lines in Seyfert 
1 type nuclei due to the reprocessing of the central hard 
radiation (CoUin-Souffrin 1987). Evidence for a disk repro- 
cessing is given by reverberation mapping studies, which 
show a very small BLR radius (several light days), in a 
contradiction with the standard cloud models predictions 
(Osterbrock 1993 and references therein). Detailed anal- 
ysis by Collin-Souffrin & Dumont (1989) shows that the 
low ionization lines (Balmer lines, for instance) from the 
disk could be successfully modeled, as this emission de- 
pends on the amount of hard X-ray radiation intercepted, 
instead of on the ionization parameter which is gener- 
ally unknown. Two types of illuminating source geometry 
have been proposed (Dumont & Colhn-Souffrin 1990a): a 
point source, located above/below the disk plane and a 
diffuse hot medium, reflecting hard radiation to the disk 
surface. The line profiles from the disk are broad, with 
FWHM ~ 10'^~'*kms~^, double-peaked in most cases, 
symmetric and non frequency shifted (Dumont & Collin- 
Souffrin 1990b). They depend mostly on disk structure, 
geometry and power of the ionizing source, as well as on 
the point of view. 

All these profile calculations have been performed un- 
der a default assumption that the disk is a planar, ax- 
isymmetric structure, illuminated uniformly. This, how- 
ever, might not always be true. If the accreting matter 
falls onto a rapidly rotating (Kerr) black hole and the an- 
gular momenta of the hole, and of the accreting gas, are 
not aligned, the disk formed is a nonplanar structure (a 
twisted or warped disk). This is the well known Bardeen- 
Petterson effect (Bardeen & Petterson 1975, Macdonald 
et al. 1986). This effect is due to the differential Lense- 

Thirring precession (wlt = — jrp™, a is the dimension- 
less spin parameter of the hole, < a < 1) of orbits around 
the spin axis of a rotating black hole (Misner et al. 1973). 
Although wlt depends on the position within the disk, 
the disk is a steady structure as the viscosity prevents 
its disassembling. Near the hole, the flow is aligned with 
the equatorial plane of the black hole, while at larger dis- 
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tanccs it is smoothly tilted to its initial orbital plane. The 
alignment radius (the Bardeen-Petterson radius) was orig- 
inally been assumed to be of the same order as the radial 
distance, where the precession time scale (~ '^'lt) com- 
parable to the infall time scale (~ RV^'^, Vr is the radial 
velocity of the gas). When the internal hydrodynamics of 
the disk is fully taken into account it appears to be much 
smaller (Kumar & Pringle 1985). 

Accretion disks spin up nonrotating, or slowly rotating, 
black holes becaiisc the angular momentum per mass unit 
of the accreting gas at the innermost stable orbit exceeds 
that of the hole (Misner et al. 1973; Shapiro & Teukolsky 
1983). The Blandford-Znajek process (electromagnetic ex- 
traction of the black hole spin energy) is probably not ef- 
ficient enough to reduce completely this accumulated spin 
momentum (Modersky et al. 1997, Livio et al. 1998). In 
other words, an AGN black hole is most probably fast ro- 
tating (o ~ 1). As a result, the accretion disk in AGN 
should be nonplanar in most cases, as there is no reason 
to suppose that the spin momentum of the accreting mass 
and the black hole spin momentum are always aligned. 
Evidence for the existence of nonplanar disks in AGN has 
been recently found by Nishiura et al. (1998). Obviously, 
such a nonplanar geometry of the disk makes it possible 
for the central radiation to reach the outer parts, as these 
can be directly "seen" from the center (Petterson 1977). 
The covering factor of a nonplanar disk in case of central 
irradiation is close to its initial inclination angle, measured 
in An units. 

In this paper, the broad Balmer emission due to repro- 
cessing of the central high energy radiation by a warped 
accretion disk is investigated. Here we present H (3 profiles. 
The profiles of other strong low ionization lines (Ha, for 
instance) could be slightly different because of different 
reprocessing properties of the medium. Our main purpose 
here is to demonstrate the effect of disk twisting on the 
line profiles, without making fits to observational data. 
That is why we choose the simplest model of a disk with 
small inclination and a point X-ray source located near the 
center. In the next sections we describe our method (Sect. 
2) and results (Sect. 3). Discussion, comparison with the 
observations, and conclusions are given in Sect. 4 

2. Modeling line profiles 

In our model a viscous, geometrically thin, and warped - 
because of the Bardeen-Petterson effect - accretion disk is 
irradiated by a point-like central source. We use a cylin- 
drical coordinate system, centered at the black hole (Fig. 
1). The irradiating source is located at Zs = 10 Rq above 
the disk plane, along the black hole axis instead of at 
the exact center, as this is a more realistic situation. The 
profiles of H/9 emission as a result of the reprocessing of 
the central hard X-ray radiation are obtained. The disk 
thickness is neglected in the computations - it does not 
affect the profile shapes significantly as the semithickness 




Fig. 1. The twisted accretion disk a steady system of 
inclined rings, each defined by two Eulerian angles and the 
radial distance (R). The axis of the black hole is aligned 
with the Z direction. To explore the irradiation, the disk is 
divided into ~ 10^ elements. Cylindrical coordinates are 
used to express the position of each element. 

to radius ratio is usually much less than the tilt angle. All 
relativistic corrections are also neglected, because they are 
small while our goal is to obtain qualitative results only. 
The disk is assiimcd to be optically thick to visual light 
and X-rays - no emission from the lower surface (probably 
also illuminated) can reach an observer located above the 
X-Y plane. 

2.1. Geometry of nonplanar disks 

Since the precession velocity is small compared to the Ke- 
plerian orbital velocity Vk = cRq^R~'^'^, the disk is a sta- 
tionary structure and can be treated as being composed 
of concentric rings, laying in different planes. Warp waves 
can propagate through the disk surface if a << 1 (Pa- 
paloizou & Pringle 1983), where a is the dimcnsionless 
viscosity parameter (Shakura & Sunyaev 1973), but this 
is not the case for AGN accretion disks, where a is usu- 
ally assumed to be 0.1-1. Each disk ring is defined by two 
Eulerian angles (3 and 7 (Fig. 1) and its radius R. 

For a stationary twisted accretion disk, f3 and 7 are 
slowly varying functions of the radial distance R: = 
/3{R), 7 = 7(i?) and p\n^oc Po, 0\r^o 0, 

7|r^oo — > 0, as is originally shown by Bardeen & Pet- 
terson (1975). In this paper we use the steady solution of 
Scheuer & Feiler (1996), derived following Pringle (1992). 
This solution, shown in Fig. 2, can be presented analyti- 
cally by: 

7(i?) « Fi?-™ 

p{R) ^ /?oe-^(«). (1) 
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Fig. 2. Eulerian angles /3 (solid line) and 7 (dashed line) 
as functions of radial distance, according Scheuer & Feiler 
(1996). Here F = 130. /3 is normalized to the initial tilt 
angle - Pq. 7 is measured in radians. 

Here the dimensionless parameter F ~ 400 -^/oa if the 
scmithickncss-to-radius ratio of the accretion disk is about 
0.01 (Collin-Souffrin & Dumont 1990), m = 0.75, and R is 
measured in Rq. A similar solution has been obtained by 
Bardccn & Pcttcrson (1975) and Hatchctt at al. (1981), 
althoughthey did not take into account the internal hydro- 
dynamics of the flow. In that case though, F and m are 
different and depend mostly on the radial velocity and 
the Kerr parameter (a), and the resulting alignment ra- 
dius is much larger. Eqs. 1 are valid for a ~ 0.3 (Kumar 
& Pringlc 1985, Natarajan & Pringlc 1998), large radial 
distances and small initial inclination angles of the flow - 
Po < 15°. If the initial tilt is significant, the influence of 
a significant amount of intercepted ionizing continuum on 
the disk structure should be taken into account and the 
disk shape should be studied numerically. 

2.2. Illumination and line profiles 

In order to explore the nonplanar disk illumination and 
the resulting line profiles we created a simple code. The 
input parameters are: the initial inclination angle /3o, the 
angles defining the line of sight, the ionizing luminosity 
Lx, the source height Zg and the black hole mass Mbh- 
For small inclinations of the disk, the deviation angle 
h = /3(i?)sin(-0 — 7(-R)) of a disk point (P) is also small 
< P < Po < 0.15) and the following approximation 
is valid: sin(/i) ~ h; cos(/i) ~ 1; see Fig. 1. The disk is 
divided into about lO'' elements, each defined by its di- 
mensions dR, dijj and coordinates R, tp. An element of the 



grid constructed in this way absorbs an amount of the cen- 
tral hard radiation proportional to its cross section area 
dSp{R,ip) (Fig. 1), which is roughly: 



d5p = dV;dii 



(l + (^ 



For all elements the following value is calculated: 

dLufsiR, V) = 4.iR^+zi) fMR, i>) 



(2) 



In the calculation of the total H P luminosity, shadow- 
ing of the elements by the inner parts of the disk must 
be taken into account. This has been done as follows: 
The contribution of each clement to the total H P lu- 
minosity from the disk is equal to dLn/S (Eq- [2]) only 
in case that there does not exist another element with 
(/ii(i?i < R,ip)- Zs/Ri) > {h{R,ijj) - Zs/R) and is in 
the opposite case. Lx - the hard X-ray luminosity, is close 
to the bolometric luminosity if the continuum extends up 
to 100 kcV with a spectral index close to 1. /H/3(i?,'0) is 
the line response to the ionizing flux, i.e. it represents the 
probability that an absorbed X-photon produces an H /?- 
photon. It depends on the physical parameters of the emit- 
ting layers, taking into account the saturation of the lines 
(due to the limited number of line-emitting atoms), the 
absorbed fraction of hard radiation (if the column density 
of the absorbing gas N < lO^^cm"^, not all hard radia- 
tion is absorbed), etc. It also may depend slightly on the 
incident flux - the reprocessing is not necessarily a lin- 
ear process (Eq. [2]). In our computations we used this 
dependence as given by Collin-Souffrin & Dumont (1989, 
1990b). In their work f-H_p{R,ijj) is presented in a multi- 
functional form, taking into account all effects mentioned. 

To obtain the profile of the line emission from the 
whole disk we have to calculate the Dopplcr shift of each 
element, i. e. the projection (Vd) of the Keplerian velocity 
onto the line of sight. Vb is given approximately by: 

Vd = VK{sm{ip - ■)/'o)cos(/io) -I- Pcos{ijj - ■y)sm{ho)) , 

where ho and ipo are the angles defining the direction to 
the observer. The whole profile can be constructed by 
adding contributions of all elements, taking into account 
their Doppler shifts. 



3. Results 

The main difference between the line emission from planar 

disks and warped ones is that in the latter case only a part 
of each disk ring is irradiated and respectively emits lines. 
A schemathic 2D map of the disk emission flux is given 
in Fig. 3. Profiles for different lines of sight are presented 
in Fig. 4. Here the bolometric luminosity is lO^'^ergs"^, 
the black hole mass is 10^ Mq , and the initial inclination 
angle is 10°. The Kerr parameter a = 1 and the viscosity 
parameter a = 0.1, so F « 130. 
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Fig. 3. A schematic map of the H /? - flux from an ir- 
radiated warped disk (Fig. 3a, b). The flux values are 

in relative logarithmic units. f3o = 10°. Mbh = 10** Af©, 
Lx = lO'^'^ergs-i, T = 130. Maps for two values of Zq, 
Zs = 10 Rg (Fig. 3a) and 100 Rq (Fig. 3b) are shown. The 
inner (i? < 10^ -Rg) part of the disk is presented in the fig- 
ure. The height (in Rq) of the disk surface in respect to 
the X-Y plane is shown in Fig. 3c. 



The most important result of this work is that the pro- 
files of the disk emission (Fig. 4) are nonsymmetrical and 
red or blue frequency shifted, as expected if only about a 
half of the disk emits in the direction towards the observer 
(Fig. 3). Only for special viewing angles, these profiles 
could be almost symmetric, double or single-peaked and 
non frequency-shifted. In Fig. 5 the profile dependence of 
various quantities is presented. T is an important parame- 
ter, affecting significantly the disk shape and respectively 
the line profiles. In Fig. 5a the profiles for F = 40 (for 
instance a = 0.1, a = 0.1); 130 (a = 1, a = 0.1) and 400 
(the upper limit - a = 1, a = 1) are shown. In Fig. 5b 
the profile dependence on the source height Zg is shown. 
Increasing Zs, a larger part of the central disk region is 
irradiated and the profiles approach those from the planar 
disks the second, symmetrically displaced peak increases 
its intensity. We do not see any physical reason, however, 
to put the source far above the central object. The de- 
pendence of the profiles on the hard X-ray luminosity is 
also significant, because of the effect of the saturation of 
the lines, which limits the emission line flux at small dis- 
tances (high velocities), independent of an increase of Lx. 
For higher Lx, the lines should become narrower (Fig. 
5c). One should keep in mind, however, that Lx should 
be 10''''^''''' erg s^^ to match the observed line intensities. 
Profile shapes are almost independent of the initial disk 
tilt (Fig. 5d). The FWHM and the shift slightly decrease 
with increasing tilt. Note that the profiles are also not de- 
pendent directly on the black hole mass and the accretion 
rate (of course, the mass and the accretion rate may af- 
fect the incident fiux at a given distance and fu_i3{R,'^) 
respectively) . 

Weaker emission from the opposite side of the disk 

may appear in case the disk is not fully opaque at visual 
wavelengths. This might be the case for the outer regions, 
where the optical depth is probably not significant if no 
dust is present (CoUin-Souffrin & Dumont 1990). Profiles 
should then be double-peaked and more or less symmetric. 
Although this case is not considered here, one can easily 
reproduce such profiles. A similar situation occurs when 
the disk is transparent to the hard radiation, which is then 
absorbed within the whole vertical structure of the disk. 

4. Discussion and conclusion 

The hard X-ray radiation, coming from the center, can 
illuminate the outer parts of the nonplanar disk, where 
optical emission lines will be emitted as a result of repro- 
cessing. Profiles of such lines arc modeled here (Fig. 4). 
Changing some basic parameters, such as the hard X-ray 
luminosity Lx, the source height .^S) F, etc., a large va- 
riety of asymmetric profiles, displaced in frequency from 
the systemic velocity profiles can be produced. 

Signatures of warped disk profiles can be seen in many 
objects (for instance 3C 227, Mkn 668, 3C 390.3, etc., Er- 
acleous & Halpern 1994). It is possible that these objects 
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Fig. 4. Line profiles of H emission from the disk, shown 
in fig. 3a, for different viewing angles. The observer is 
tilted at an angle ho = 45° with respect to the black hole 
equatorial plane. The azimuthal angle along the X-Y plane 
{tpo) is respectively 0,30, ..150° - fig. 4a. .f. The inner and 
the outer radii of the disk are 10 and 10* Rq respectively. 
Abscissa is in lO^kms"^. 



contain an illuminated warped disk, but of course, other 
explanations can not be excluded. However, most of the 
AGN show symmetric profiles, without frequency shifts. 
On average, no more than roughly 10-15% of the broad 
line emission in AGN can be reproduced by an irradiated 
nonplanar disk. Wc note that this value is close to the cov- 
ering factor of such a disk. One may conclude, therefore, 
that the bulk of the line emission in AGN arises from 
some more or less spherical or conical structure (a sys- 
tem of clouds or star atmospheres, a jet), with a covering 
factor which is close to unity, but not from a irradiated 
thin disk. In that case, the presence of an inclined disk 
structure will cause only asymmetry of lines or additional 
peaks, displaced in frequency from the main profile. There 
are several possibilities for objects where no such asym- 
metries or displaced peaks are observed. 
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Fig. 5. The profile dependence for different factors: the 
value of r fig. 5a; the source height fig 5b (Zg is 
measured in Rq)] the central hard luminosity - fig. 5c 
(log(Lx): with Lx in ergs^^) and the initial tilt angle - 
fig. 5d (/9o is in degrees), ho — 45°, ^po = 0° . Abscissa is 
in lO^kms"^. All profiles are normalized to the maximal 
value. 

i) The presence of a thin disk in AGN is quite unusual. 
The disk may also be thick (slim, advection dominated 

Rces 1984; Blandford & Begelman 1998; Chakrabarti 
1998; Park & Ostriker 1998). 

ii) The central black hole is nonrotating and the disk is 
not warped. In this case, if the irradiation is possible, the 
profiles should be symmetric and double-peaked, which 
is also seldom observed, however. Similar profiles should 
be observed in case of a disk transparent to visual light 
and/or X-rays. 

iii) The spins of the infalling gas and the hole are 
aligned and the disk is not warped. This is possible if 
the accreting matter is supplied from a single direction 
for a long period of time. As a result the black hole would 
align its spin with the spin of the accreting matter on rel- 
atively short time scale - 10^~^ yr (Scheuer & Feiler 1996; 
Natarajan & Pringle 1998). 

iv) The disk is warped, but the illuminating source 
is positioned far above the disk plane - at ^ 10^ Rq or 
more (or a huge reflecting corona with a significant optical 
depth is present). In this case a small inclination of the 
disk would not affect significantly the line profiles. 

Nevertheless, many objects reveal asymmetries and 
displaced peaks, which could be successfully reproduced 
by a warped accretion disk. These profiles depend strongly 
on the twisting structure, the illuminating source geom- 
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etry and the line of sight to the observer. In principle 
they can be modeled with high signal-to-noise spectra of 
the object by varying the model parameters. Knowledge 
of the twisting structure may allow the determination of 
some very important characteristics of AGN, at first or- 
der the spin momentum of the black hole and the viscosity 
parameter of the accretion disk. Their determination by 
other methods is still unreliable. 
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